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1. Introduction 

After the discovery of a Higgs boson at the LHC 
in 2012 (HE], one of the most important remaining 
questions is whether the newly discovered particle 
is part of an extended scalar sector. A CP-odd 
Higgs boson, A, appears in many models with an 
extended scalar sector, e.g. in the case of the two- 
Higgs-doublet model (2HDM) [3]. 

The addition of a second Higgs doublet leads to 
five Higgs bosons after the electroweak symmetry 
breaking. The phenomenology of such a model is 
very rich and depends on the vacuum expectation 
values of the Higgs doublets, the CP properties of 
the Higgs potential and the values of its parame¬ 
ters and the Yukawa couplings of the Higgs dou¬ 
blets with the fermions. In general, it is possible 
to accommodate in the model a Higgs boson com¬ 
patible to the one discovered at the LHC. In the 
case where the Higgs potential of the 2HDM is CP- 
conserving, the Higgs bosons after electroweak sym¬ 
metry breaking are two CP-even {h and H), one 
CP-odd ( A ) and two charged (H ± ) Higgs bosons. 
Many theories beyond the Standard Model (SM) in¬ 
clude a second Higgs doublet, such as the minimal 
supersymmetric SM (MSSM) dH9, axion models 
(e.g. Ref. w and baryogenesis models (e.g. Ref. 

ITS! j. Searches for a CP-odd Higgs boson are re¬ 


ported in Refs. mm- 

In this Letter, a search for a heavy CP-odd Higgs 
boson decaying into a Z boson and the ~125 GeV 
Higgs boson, h, is described. The A -A Zh decay 
rate can be dominant for part of the 2HDM pa¬ 
rameter space, especially for an A boson mass, ipa, 
below the tt threshold. In this case, the A boson 
is produced mainly via gluon fusion and its natural 
width is typically small: Pa/uia < 0(1%). 

The search is performed for niA in the range 220 
to 1000 GeV, reconstructing Z -» ££ decays (where 
£ = e, \x) with h —► bb or h —> tt , as well as Z —► vv 
with h —► bb. The selected h boson decay modes 
provide high branching ratios and the possibility to 
fully reconstruct the Higgs boson decay kinemat¬ 
ics. The reconstructed invariant mass (or trans¬ 
verse mass) of the Zh pair, employing the measured 
value of the h boson mass, to^, to improve its res¬ 
olution, is used to search for a signal. 

2. Data and simulated samples 

The data used in this search were recorded with 
the ATLAS detector in proton proton collisions at 


1 Throughout this Letter, the notation h —> bb, h —» tt, 
Z —^ w and Z —» tl is used for h —> bb, h —> t + t~, Z —» vv 
and Z —> (A i~, respectively. 
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a centre-of-mass energy of 8 TeV. The ATLAS de¬ 
tector is described in detail elsewhere JT5]. The 
integrated luminosity of the data sample, selecting 
only periods where all relevant detector subsystems 
were operational, is 20.3 ± 0.6 fb -1 [TB]. The data 
used in the IAt and Ubb final states were collected 
using a combination of single-electron, single-muon, 
dielectron (ee) and dimuon (p/i) triggers. Depend¬ 
ing on the trigger choice, the p^\ thresholds vary 
from 24 to 60 GeV for the single-electron and single¬ 
muon triggers, and from 12 to 13 GeV for the ee and 
pp triggers. The data used in the vvbb final state 
were collected with a missing transverse momentum 
(E^ iss ) trigger with a threshold of E^ iss > 80 GeV. 

Signal events from a narrow-width A boson 
produced via gluon fusion are generated with 
MadGraph5 [T7j for all final states considered in this 
search. The parton showering is performed with 
PYTHIA8 PUGH. 

Production of W and Z bosons in association 
with jets is simulated with SHERPA [20]. Top-quark 
pair and single top-quark production is simulated 
with POWHEG I5TH53] and AcerMC [21]. Production 
of WW, WZ, and ZZ dibosons are simulated us¬ 
ing POWHEG. The WZ and ZZ processes include the 
production of off-shell Z bosons (Z*) and photons 
( 7 *). Triboson production (WWW«, ZWW<*\ 
ZZZ (*)) and top pair production in association 
with a Z boson are generated with MadGraph5. 
Finally, the production of the SM Higgs boson 
in association with a Z boson is considered as a 
background in this search. It is simulated using 
PYTHIA8. 

The CTEQ6L1 [25] set of parton distribution 
functions was used for samples generated with 
MadGraph5 and PYTHIA 8 . The CT10 [26] set was 
used for the other samples. 

All generated samples are passed through the 
GEANT4-based [2?( detector simulation of the AT¬ 
LAS detector [2%] . The simulated events are over¬ 
laid with minimum-bias events, to account for the 
effect of multiple interactions occurring in the same 
and neighbouring bunch crossings (“pile-up”). The 


2 ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point (IP) in the centre of 
the detector and the z-axis along the beam pipe. The x-axis 
points from the IP to the centre of the LHC ring, and the 
y -axis points upward. Cylindrical coordinates (r, 0) are used 
in the transverse plane, <j> being the azimuthal angle around 
the beam pipe. The pseudorapidity is defined in terms of 
the polar angle 9 as q = — lntan(0/2). Transverse momenta 
are computed from the three-momenta, p, as p-y = \p] sin 9. 


events are reweighted so that the average number 
of interactions per bunch crossing agrees with the 
data. 

The background estimation in this search for 
most processes is based on data driven techniques, 
but in some cases only simulated samples are used. 
In that case, the simulated samples are normalized 
using theoretical cross section calculations. In par¬ 
ticular, for diboson production both qq [29] and 
gg PSE .31. initiated processes are included. Tri¬ 
boson production follows Ref. [35] and top pair 
production in association with a Z boson follows 
Refs. [331 151] . SM Higgs boson production in asso¬ 
ciation with a Z boson uses a calculation described 
in Ref. [85] . 

3. Object reconstruction 

Electrons are identified from energy clusters in 
the electromagnetic calorimeter that are matched 
to tracks in the inner detector [36]. Electrons are re¬ 
quired to have \rj\ < 2.47 and pt > 7 GeV. Isolation 
requirements, defined in terms of the calorimetric 
energy or the px of tracks within cones around the 
object, as well as quality requirements are applied 
to distinguish electrons from jets. 

Muons are reconstructed by matching tracks re¬ 
constructed in the inner detector to tracks or track 
segments in the muon spectrometer systems [37 . 
The muon acceptance is extended to the region 
2.5 < \p\ < 2.7, which is outside the inner detec¬ 
tor coverage, using only tracks reconstructed in the 
forward part of the muon detector. Muons used for 
this search must have |p| < 2.7, px > 6 GeV and 
are also required to pass isolation requirements. 

Jets are reconstructed using the anti-Ay algo¬ 
rithm [38] with radius parameter R = 0.4 and px > 
20 GeV (p T > 30 GeV) for |p| < 2.5 (2.5 < 
H < 4.5). Low-px jets from pile-up are rejected 
with a requirement on the scalar sum of the px 
of the tracks associated with the jet: for jets with 
I 77 I < 2.4 and px < 50 GeV, tracks associated with 
the primary verte 30 must contribute over 50% to 
the sum. Jets from the decay of long-lived heavy- 
flavour hadrons are selected using a multivariate 
tagging algorithm ( 6 -tagging) [35]. The 6 -tagging 
efficiency is 70% for jets from 6 -quarks in a sample 
of simulated tt events. 


2 The primary vertex is taken to be the reconstructed ver¬ 
tex with the highest Sp;], of the associated tracks. 
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Hadronic decays of r leptons (rhad) 00] are re¬ 
constructed starting from clusters of energy in the 
calorimeter. A Thad candidate must lie within 
|7/| < 2.47, have a transverse momentum greater 
than 20 GeV, one or three associated tracks and 
a total charge of ±1. Information on the collima- 
tion, isolation, and shower profile is combined into 
a multivariate discriminant to reduce backgrounds 
from quark- or gluon-initiated jets. Dedicated al¬ 
gorithms that reduce the number of electrons and 
muons misidentified as hadronic r decays are ap¬ 
plied. In this analysis, two Thad identification se¬ 
lections are used —“loose” and “medium”— with 
efficiencies of about 65% and 55%, respectively. 

The missing transverse momentum (_Kj? lss ) is 
computed using fully calibrated and reconstructed 
physics objects, as well as clusters of calorimeter¬ 
cell energy deposits that are not associated with 
any object 02- In addition, a track-based missing 
transverse momentum ( p is calculated as the 
negative vector sum of the transverse momenta of 
tracks with |p| < 2.4 and associated with the pri¬ 
mary vertex. 


4. Search for A —> Zh with h —> tt 

In the search for A —> Zh —» ttrr, three chan¬ 
nels are considered, distinguished by the way the rr 
pair decays: two r leptons decaying hadronically 
(ThadThad), one leptonic and one hadronic decay 
(TiepThad) and, finally, two leptonic decays (T lep Ti e p). 
Electrons in the Thad^had and Ti ep Th a d channels are 
rejected in the transition region between the barrel 
and end-cap of the detector (1.37 < |7y| < 1.52). 
Muons in the Th a dTh a d and Ti ep Th a d channels are 
considered only for |?y| < 2.5. 

The resolution of the reconstructed A boson mass 
is improved using a mass-difference variable, 

m A c = rn UTT - mu - rn TT + m z + m h , 

where m z is the mass of the Z boson, mh = 
125 GeV is the mass of the CP-even Higgs bo¬ 
son, mu is the invariant mass of the two leptons 
associated with the Z boson decay, and m« rr de¬ 
notes the Htt invariant mass. The value of m TT , 
the invariant mass of the r’s, is estimated with the 
Missing Mass Calculator (MMC) )42| . The mass 
resolution for all rr channels ranges from 3% at 
m A = 220 GeV to 5% at m A = 1 TeV. 


4 -l ■ tip ii H d ri,ad 

Events in the tirhadThad channel are required to 
contain exactly two opposite-sign leptons It (ee or 
pp) and exactly two opposite-sign Th a d- The px 
requirements for these objects are px > 26 GeV 
(15 GeV) for the leading (subleading) electron, 
Pt > 25-36 GeV (10 GeV) for the leading (sub¬ 
leading) muon, depending on the trigger, and pt > 
35 GeV (20 GeV) for the leading (subleading) rhad 
candidates. The rh a d candidates are required to 
satisfy the “loose” rhad identification criterion. In 
addition, the ee/pp invariant mass and the rr in¬ 
variant mass have to lie in the ranges 80 < mu < 
100 GeV and 75 < m TT < 175 GeV. Finally, the pt 
of the it pair, p^, is required to be: 

125 GeV, if mj c > 400 GeV 
0.64 x m] 4 C — 131 GeV, otherwise. 

This requirement maximizes the sensitivity over the 
whole explored A mass range. In the region of 
pf > 125 GeV, there is little background present, 
so tightening the requirement results in no ad¬ 
ditional increase in sensitivity. The total accep¬ 
tance times selection efficiency varies from 6.2%, 
for m A = 220 GeV, to around 18% for the highest 
A boson masses considered. 

The dominant background for this channel origi¬ 
nates from events where one or both of the Thad’s is 
a misidentified jet (“fake-Th a d background”). This 
background is dominated by Z+ jets events, with 
small contributions from dibosons and events with 
top quarks, and it is estimated using a template 
method. The shape of the fake-rh a d background is 
taken from a control region (the “template region”) 
that contains events satisfying all the tirhadThad se¬ 
lection criteria apart from the requirements for an 
opposite-sign Th a d T had pah' and the Thad identifica¬ 
tion criteria. The fake-Thad background is normal¬ 
ized by using two additional control regions. The 
first region, “A”, contains events that satisfy the 
signal selection criteria, with the exception that the 
m TT constraint is inverted, i.e. m TT < 75 GeV or 
m TT > 175 GeV. The second region, “B”, contains 
events that satisfy all the template selection crite¬ 
ria, with the exception that the m TT constraint is 
inverted, as in the region “A” definition. The ratio 
of the number of events in “A” to the number of 
events in “B” is used to scale the template region 
events in order to obtain the normalization of the 
fake-Th a d background. 

In addition to the fake-Thad background, there 
are also contributions from backgrounds with real 
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•^Thad'Hiad objects in the event. These backgrounds 
come primarily from ZZ^ productiorQ SM Higgs 
boson production in association with a Z boson is 
estimated using simulation, and contributes 17% of 
the total background. 

J^.2. ^TlepThad 

Events in the £tTi ep Th a d channel are required to 
contain exactly three light leptons, ppp, epp, eep 
or eee, and exactly one Thad- The px requirements 
for these objects are px > 26 GeV (15 GeV) for the 
leading (remaining) electron(s), px > 25-36 GeV 
(10 GeV) for the leading (remaining) muon(s), de¬ 
pending on the trigger, and px > 20 GeV for the 
Ti ia d ■ Subsequently, all the possible U pairs that are 
composed of opposite-sign, same-flavour leptons are 
selected. From these pairs, the pair that has the 
invariant mass closest to mz is considered to be 
the lepton pair from the Z boson decay. The third 
light lepton is considered to be the leptonic r de¬ 
cay, and it is used along with the Thad to define the 
PepThad pair. This light lepton is required to have 
opposite-sign charge with respect to the Thad- Hr ad¬ 
dition, the Thad is required to satisfy the “medium” 
Th a d identification requirement, and mu and m TT 
have to lie in the ranges 80 < mu < 100 GeV 
and 75 < m TT < 175 GeV. The total accep¬ 
tance times selection efficiency varies from 6% for 
toa = 220 GeV, to around 17% for the highest A 
boson masses considered. 

About half of the total background for this chan¬ 
nel comes from events where the Thad and/or the 
light lepton is a misidentified jet (“fake-r/i back¬ 
ground”). This background is dominated by di¬ 
boson and Z +jets events and it is estimated us¬ 
ing a template method. The shape of the fake- 
T/t background is taken from a control region (the 
“template region”) that contains events satisfy¬ 
ing all ^fViep T had selection criteria, apart from re¬ 
quiring “medium” Thad identification criterion and 
opposite-sign charge for the Ti ep Th a d pair. The fake- 
T/i background is normalized by using two addi¬ 
tional control regions, defined similarly to those in 
the ^ThadThad channel. 

The other half of the background comes from 
events with real ^TiepThad objects in the event. 
These backgrounds come primarily from ZZ^*4 pro¬ 
duction. There is also a small (11%) contribution 


4 The notation ZZZ'Z is used here to include ZZ, ZZ* 
and Z 7 * 


from the SM Higgs boson production in association 
with a Z boson, which is estimated using simula¬ 
tion. 

4 ■ A f f T| ep T]ep 

Events in the ilT\ op T\ ep channel are required to 
contain at least four leptons, which form one same- 
flavour and opposite-sign pair consistent with the 
Z mass (80 < mu <100 GeV), and either a same- 
flavour or different-flavour pair with an invariant 
mass reconstructed with the MMC algorithm, con¬ 
sistent with a decay from the CP-even Higgs bo¬ 
son ( 90 < m TT < 190 GeV). One muon is al¬ 
lowed to be reconstructed in the forward region 
(2.5 < |p| < 2.7) of the muon spectrometer, or to 
be identified in the calorimeter with px > 15 GeV 
and 1 77 1 < 0.1 [37]. The highest-px lepton must sat¬ 
isfy px > 20 GeV, and the second (third) lepton in 
px order must satisfy px > 15 GeV (px > 10 GeV). 
Among all the possible lepton quadruplets in an 
event the one minimizing the sum of the mass dif¬ 
ferences with respect to both the Z and h bosons 
is chosen. 

Two different analysis categories are defined 
based on the lepton flavours in the Higgs boson de¬ 
cay: ee or pp (SF), and ep (DF). The expected 
background is very different in the two cases. For 
the SF channel, the background is dominated by 
ZZ M production with Z ee/pp decays. For the 
DF channel, the main background is from the ZZ W 
process through the Z — > Ti op Ti op decay chain, but 
other backgrounds are also important. The signal- 
to-noise ratio in the SF category is improved by 
using a set of requirements specifically targeted to 
suppress the main ZZ background. First, a veto 
on the on-shell production of Z boson pairs is intro¬ 
duced, requiring the invariant mass of the h boson 
leptons to lie outside the Z peak: < 80 GeV 

or mh > 100 GeV. Background events are charac¬ 
terized by low missing transverse momentum and 
are further rejected by requiring E™ lss > 30 GeV, 
and the azimuthal angle between the E“ lss direc¬ 
tion and the Z boson transverse momentum to be 
greater than tt/2. Furthermore, a requirement that 
the highest-px lepton of the U pair associated with 
the h boson has px > 15 GeV is applied, since it 
is found to be effective against backgrounds from 
Z+ jets production. The total acceptance times se¬ 
lection efficiency varies from 6.5% (1.5%) for DF 
(SF) channel for = 220 GeV, to around 20% 
for both channels for the highest A boson masses 
considered. 
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The subleading contributions to the background 
are from diboson and triboson production, tt pro¬ 
duction in association with a Z boson, and SM 
Higgs boson production. All these are determined 
from simulation and amount to about 95% (65%) 
of the total background in the SF (DF) category. 
The other background events have at least one lep¬ 
ton which is a misidentified jet or a lepton from 
a heavy-flavour quark decay and are dominated 
by Z+jets production, with a smaller contribution 
from top-quark production. These backgrounds are 
estimated using a control region where one or both 
of the leptons in the it pair associated with the 
h —> riepTiep decay fail to satisfy the isolation cri¬ 
teria. After subtraction of genuine sources of four- 
lepton events using simulation, the data are extrap¬ 
olated to the isolated signal region using normaliza¬ 
tion factors derived from simulated samples. 

4-4- Systematic uncertainties and results 

The most important systematic uncertainty for 
the backgrounds with real Utt objects in the 
TiepThad and Ti ep Ti op channels comes from the un¬ 
certainty on the theoretical cross sections used in 
the normalization. They are due to the parton dis¬ 
tribution function choice, the renormalization and 
factorization scales, as well as the as value. This 
amounts to an uncertainty on the normalization of 
this background of about 5.0% for the TiepThad chan¬ 
nel and 6.4% for Ti ep Ti op . In the ThadThad channel, 
the largest contributions come from the Th a d identi¬ 
fication and energy scale and amounts to 8.9% [30] , 
The fake- Thad background systematic uncertainty 
for the tt channels is dominated by the statisti¬ 
cal uncertainty on data in control regions used for 
the background normalization. It amounts to a 
normalization uncertainty of 38% and 25% for the 
TiepThad and ThadThad channels, respectively. For 
the Ti e pTi e p channel, the normalization uncertainty 
is 65% (25%) for the SF (DF) category. 

The reconstructed A boson mass distributions 
for events passing the ^ThadThad, ^TiepThad and 

Ti e pTiep selections are shown in Fig. [l] The num¬ 
ber of events passing the Utt channel selections are 
shown in Table[l] The agreement of the expectation 
with data is very good. 

5. Search for A —* Zh with h —y bb 

This section describes the searches in the A —► 
Zh —> iibb and A —>• Zh —> vvbb channels. 


Table 1 

The number of predicted and observed events for 
the Utt channels. 



Expected Background 

Data 

^'ThadThad 

28 ±6 

29 

^^"lepThad 

17 ±4 

18 

^TiepTiep (SF) 

9.5 ±0.6 

10 

^T^Tiep (DF) 

7.2 ±0.7 
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5.1. iibb selection 

Events in the iibb channel are selected by requir¬ 
ing either two electrons or two muons. In the case 
of muons they are required to be of opposite-sign 
charge. Leptons must have px > 7 GeV, and elec¬ 
trons are restricted to \rj\ < 2.47, while muons must 
have | r ]| < 2.7. Tighter acceptance requirements are 
placed on one of the leptons in each event in order 
to select a sample for which the trigger efficiency is 
high and to reduce the nrulti-jet background, while 
keeping a high signal acceptance. These require¬ 
ments are that the leptons have pt > 25 GeV, and, 
if they are muons, satisfy |p| < 2.5. A dilepton 
invariant mass window of 83 < mu < 99 GeV is 
imposed to reduce top-quark and multi-jet back¬ 
grounds. 

The h —» bb decay is reconstructed by requiring 
two 6-tagged jets with px > 45 GeV (20 GeV) for 
the leading (subleading) jet. Events with more than 
two 6-tagged jets are removed but all events with 
one or more additional jets failing 6-tagging are re¬ 
tained. The h —>• bb decay is selected by requiring 
that the invariant mass of the two 6-tagged jets lies 
within the range 105 < nibb < 145 GeV. 

The top-quark background, which includes top- 
quark pair and single top-quark production, is 
reduced by requiring E™ lss /^/H t < 3.5 GeV 1 / 2 , 
where Ht is defined as the scalar sum of the px 
of all jets and leptons in the event. 

The reconstructed A boson mass, to(| c , is the in¬ 
variant mass of the two leptons and two 6-tagged 
jets. In this calculation, the four-momentum of 
each 6-tagged jet is scaled by 125 GeV/mbb in or¬ 
der to improve the resolution. The resulting m r J c 
resolution ranges from 2% at mA = 220 GeV to 3% 
at mA = 1 TeV. 

In order to reduce the dominant Z +jets back¬ 
ground, a requirement is imposed on the transverse 
momentum of the Z boson, p;^, reconstructed from 
the two leptons: pf > 0.44 x m —106 GeV, where 
mA is in units of GeV. The requirement depends on 


5 








Fig. 1. Distributions of the reconstructed A boson mass for the combined iiT\ VM \T\ VM \ and f^Ti ep Thad final 
states (a) and the UT\ ep T\ ep final states (b). The signal shown in both cases corresponds to a(gg — > A) x 
BR(A — > Zh) x BR(/i — > tt) = 50 fb with toa = 340 GeV. The background contributions shown are 
the results of simulation and data-driven estimation methods. The background uncertainty is shown as a 
hatched area, and the overflow is included in the last bin. 


m] J c since the background is generally produced at 
low whereas the mean pip increases with tua 
for the signal. The total acceptance times selection 
efficiency varies from 7%, for niA = 220 GeV, to 
around 16% for the highest A boson masses consid¬ 
ered. 

5.2. vvbb selection 

The event selection in the vvbb channel follows 
closely the SM h —» bb analysis in Ref. [43] . Events 
are selected with E“ iss > 120 GeV, p^ iss > 30 GeV 
and no electrons or muons with px > 7 GeV. In 
addition to the jet selection of the ££bb analysis, 
additional restrictions are applied. In order to sup¬ 
press top-quark background, which is larger than 
in the ££bb channel, events are rejected if any of 
the following conditions is satisfied: there is a jet 
with |? 7 | > 2.5; there are four or more jets; one 
of the 6 -tagged jets is the third-highest-px jet. In 
order to select a sample for which the trigger effi¬ 
ciency is high, Ht is required to be above 120 GeV 
(150 GeV) for events with two (three) jets. There 
are also requirements on the separation between 
the two 6 -jets in the rj-cj) space, A Rjj, to suppress 
Z +jets and VF+jets backgrounds as described in 
Ref. [43]. As in the ££bb channel, the h boson is 
selected by requiring 105 < m bb < 145 GeV. 

Additional requirements are imposed on angu¬ 
lar quantities sensitive to the presence of neu¬ 


trinos in order to suppress the multi-jet back¬ 
ground: the azimuthal angle between E™ lss and 
p T miss : A (f>(E™ iss ,p™ iss ) < tt/ 2; the mini¬ 

mum azimuthal angle between E™ lss and any jet 
min[A^>(.E™ lss , jet)] > 1.5; and the azimuthal angle 
between E™ lss and the 6 -jet pair A</)(E™ 1SS , 66 ) > 
2.8. The total acceptance times selection efficiency 
varies from 4%, for toa = 400 GeV, to around 7% 
for the highest A boson masses considered. 

It is not possible to accurately reconstruct the 
invariant mass of the A boson due to the pres¬ 
ence of neutrinos in the final state. Therefore, 
the transverse m ass is used as the final discrimi- 
nant: m™ c ’ T = ( E+ E ™ lss ) 2 — (p^ b + E7p lss ) 2 , 

where E^ and p^ b are the transverse energy and 
transverse momentum of the 6 -jet pair system. As 
in the ££bb channel, the resolution is improved 
by scaling each 6 -tagged jet four-momentum by 
125 GeV/m bb . 

5.3. Backgrounds 

All backgrounds in ££bb/vvbb final states are de¬ 
termined from simulation, apart from the multi¬ 
jet background, which is determined from data. 
The multi-jet background in the / 1/166 final state 
is found to be negligible. In the eebb final state, 
the background is determined by selecting a sam¬ 
ple of events with the electron isolation require¬ 
ment inverted. The sample is normalized by fitting 
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the mu distribution. In the vvbb final state, the 
multi-jet background is determined by inverting the 
A^(£’™ lss ,p 1 ? llss ) requirement. The sample is nor¬ 
malized using the region with min[A0(£J™ lss , jet)] < 
0.4. 

The Z+jets simulated sample is split into differ¬ 
ent components according to the true flavour of the 
jets, i.e. Z + 11, Z + cl, Z + cc, Z + bl, Z + bc and 
Z + bb , where l denotes a light quark (u, d, s) or a 
gluon. These components are constrained by defin¬ 
ing control samples which have the same selection 
as the Mbb final state, but with the requirements 
on the number of 6-tagged jets changed to either 
zero or one. The samples are further divided into 
events with two or at least three jets. In order to 
improve the description of the data, corrections are 
applied to the simulation as a function of the az¬ 
imuthal angle between the two leading jets, A 
for Z + ll events and a function of pf for the other 
components, as described in detail in Ref. [15] , 

The kF+jets background, which contributes sig¬ 
nificantly only in the vvbb final state, is split into its 
components in the same way as the Z+jets sample. 
It is constrained by defining a sample of events that 
are selected using the U™ lss triggers and contain ex¬ 
actly one lepton with px > 25 GeV and a tightened 
isolation requirement. The transverse momentum 
of the lepton and E™ lsa system (p^) is required 
to be above 120 GeV to approximately match the 
phase space of the signal region. The sample is split 
into events with zero, one or two 6-tagged jets and 
into events with 2 and 3 jets. A correction depend¬ 
ing on A </>jj is applied to W + 11 and W + cl events, 
following studies similar to those performed for the 
Z+jets background 051 . 

A correction is made to the px distribution of 
tt production in the simulation to account for an 
observed discrepancy with the data [44] , The nor¬ 
malization of top-quark pair production in the £tbb 
channel is measured by defining a sample of events 
with exactly one electron and one muon, one of 
which has px > 25 GeV, and two 6-tagged jets with 
50 < mbb < 180 GeV. 

5-4- Systematic uncertainties and results 

The most important experimental systematic un¬ 
certainties in the Mbb and vvbb final states come 
from the jet energy scale uncertainty and the b- 
tagging efficiency. 

The jet energy scale systematic uncertainty arises 
from several sources including uncertainties from 


the in situ calibration, pile-up dependent correc¬ 
tions and the jet flavour composition |45| . In ad¬ 
dition, an uncertainty on the jet energy resolution 
is applied. The jet energy scale and resolution un¬ 
certainties are propagated to the E™ lss . The un¬ 
certainty on E™ ss also has a contribution from 
hadronic energy that is not associated with jets m- 

The 6-tagging efficiency uncertainty depends on 
jet px and comes mainly from the uncertainty on 
the measurement of the efficiency in tt events [55] , 
Similar uncertainties are derived for the c-tagging 
and light-flavour jet tagging [36] , 

Other experimental systematic uncertainties that 
are included but have a smaller impact are uncer¬ 
tainties from lepton energy scale and identification 
efficiency, the efficiency of the E™ 1SB trigger and the 
uncertainty on the multi-jet background estimate, 
which is taken to be 100% of the estimated number 
of events. 

In addition to the experimental systematic un¬ 
certainties, modelling systematic uncertainties are 
applied, accounting for possible differences between 
the data and the simulation model used for each 
process. For the background samples, the proce¬ 
dure described in Ref. [43] is followed. The Z+jets 
and tU+jets backgrounds include uncertainties on 
the relative fraction of the different flavour com¬ 
ponents, and on the mbb, and VtIPt distri¬ 

butions. For tt production, uncertainties on the 
top-quark transverse momentum, mbb, E™ lss and 
Pt/Pt distributions are included. Uncertainties on 
the ratio of two-jet to three-jet events are also in¬ 
cluded for each background. 

The m™ c and m^ c ’ T distributions for events pass¬ 
ing the Mbb and vvbb final-state selections, respec¬ 
tively, are shown in Fig. [2] The distributions are 
shown after a profile-likelihood fit, which constrains 
simultaneously the signal yield and the background 
normalization and shape, which is performed in 
the same manner as in Ref. 03] • The overall 
background is more constrained than the individ¬ 
ual components, causing the errors of indvidual 
components to be anti-correlated. The number of 
events passing the Mbb and vvbb final state selec¬ 
tions are shown in Table [2] where the values for the 
expectations and uncertainties are obtained from 
the profile-likelihood fit. 

6. Results 

In all channels, no significant excess of events 
is observed in the data compared to the predic- 
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mj c [GeV] 



m™ J [GeV] 


(a) Ubb 


(b) uubb 


Fig. 2. Distributions of the reconstructed A boson mass for the Ubb final state (a) and the A boson 
transverse mass for the vvbb final state (b). The signal shown in both cases corresponds to cr(gg —> A) x 
BR(A —> Zh ) x BR(/i —> bb) = 500 fb with toa = 500 GeV. The predicted distributions are shown after the 
profile likelihood fit to the data. The uncertainty is shown as a hatched area, and the overflow is included 
in the last bin. 



m A [GeV] 

(a) A — > Zh , h — ¥ tt 



m A [GeV] 

(b) A —> Zh, h ->■ bb 


Fig. 3. Combined observed and expected upper limits at the 95% CL for the production cross section of a 
gluon-fusion-produced A boson times its branching ratio to Zh and branching ratio of h to (a) tt and (b) 
bb. The expected upper limits for subchannels are also shown. 


tion from SM background sources. The significance 
of local excesses is estimated using p-values calcu¬ 
lated with a test statistic based on the profile like¬ 
lihood m- The largest data excesses are at toa = 
220 GeV (p-value = 0.014) and toa = 260 GeV 
[p -value = 0.14) in the combined final states with 


h -A bb and h —► tt, respectively. Exclusion lim¬ 
its at the 95% confidence level (CL) are set on the 
production cross section times the branching ratio 
BR(A —» Zh) as a function of the A boson mass. 
The exclusion limits are calculated with a modi¬ 
fied frequentist method [45], also known as CLs, 

































(a) Type-I 2HDM, = 300 GeV 



cos(P-a) 

(c) Lepton-specific 2HDM, m a = 300 GeV 



(b) Type-II 2HDM, ttia — 300 GeV 



cos(P-a) 

(d) Flipped 2HDM, 4 = 300 GeV 


Fig. 4. The interpretation of the cross-section limits in the context of the various 2HDM types as a function 
of the parameters tan /3 and cos(/3 — a) for rriA = 300 GeV: (a) Type-I, (b) Type-II, (c) Lepton-specific, and 
(d) Flipped. Variations of the natural width up to Ta/itia=5% are taken into account. For Type-II and 
Flipped 2HDM, the 6-associated production is included in addition to the gluon fusion. The narrow regions 
with no exclusion power in Type-I and Type-II at low tan/3 and far from cos(/3 — a) = 0 are caused by 
vanishing branching ratios of h —> bb and/or h —»• rr. The blue shaded area denotes the area excluded by 
taking into account the constraints on the CP-odd Higgs boson derived by considering the A —> rr decay 
mode after reinterpreting the results in Ref. m- 9 






















m A [GeV] m A [GeV] 


(a) Type-I 2HDM, cos(/3 — a) = 0.1 


(b) Type-II 2HDM, cos(/3 — a) = 0.1 




m A [GeV] 


m A [GeV] 


(c) Lepton-specific 2HDM, cos(/3 — a) = 0.1 


(d) Flipped 2HDM, cos(/3 — a) = 0.1 


Fig. 5. The interpretation of the cross-section limits in the context of the various 2HDM types as a function 
of the parameters tan/3 and for cos(/3 — a) = 0.1: (a) Type-I (a), (b) Type-II, (c) Lepton-specific, and 
(d) Flipped. Variations of the natural width up to T a /?ti. a = 5% are taken into account. The grey solid 
area indicates that the width is larger than 5% of in a- For Type-II and Flipped 2HDM, the b- associated 
production is included in addition to the gluon fusion. The blue shaded area denotes the area excluded by 
taking into account the constraints on the CP-odd Higgs boson derived by considering the A —» rr decay 
mode after reinterpreting the results in Ref. ini- 
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Table 2 

Predicted and observed number of events for the 
££bb and vvbb final states shown after the profile 
likelihood fit to the data. 


££bb vvbb 


-Z+jets 

1443 

± 

60 

225 

± 

11 

PF+jets 


- 


55 

± 

8 

Top 

317 

± 

28 

203 

± 

15 

Diboson 

30 

± 

5 

10.8 

± 

1.6 

SM Zh, Wh 

31.7 

± 

1.8 

22.5 

± 

1.2 

Multi-jet 

20 

± 

16 

3.2 

± 

3.1 

Total background 

1843 

± 

34 

521 

± 

12 


Data 1857 511 


and the profile likelihood method, using the binned 
to^ c mass distributions for ££tt and ££bb final states 
and the binned m^ c,T distribution for the vvbb final 
state. 

Fig-! shows the 95% CL limits on the produc¬ 
tion cross section times the branching ratio, cr(gg —► 
A) x BR(A —► Zh) x BR(6 —> bb/rr ), as well as the 
expected limits for each individual subchannel. The 
limit on the production times the branching ratio 
is in the range 0.098-0.013 pb and 0.57-0.014 pb 
for via in the range 220-1000 GeV for the rr and 
bb channels, respectively. The rr channels use few 
signal mass points beyond mi = 500 GeV, since 
a coarse binning in m r J c is adopted in view of the 
very small predicted number of background events. 

The results of the search in the rr and bb channels 
are combined in the context of the CP-conserving 
2HDM [3J, which has seven free parameters and four 
arrangements of the Yukawa couplings to fermions. 
In particular, the free parameters are the Higgs 
boson masses (m?,, run, tua, m#±), the ratio of 
the vacuum expectation values of the two dou¬ 
blets (tan/3), the mixing angle between the CP- 
even Higgs bosons (a) and the potential param¬ 
eter m\ 2 that mixes the two Higgs doublets. The 
Yukawa coupling arrangements distinguish four dif¬ 
ferent 2HDM models, determining which of the two 
doublets, and 4> 2 , couples to up- and down- 
type quarks and leptons. In the Type-I model, <1> 2 
couples to all quarks and leptons, whereas in the 
Type-II, couples to down-type fermions and $ 2 
couples to up-type fermions. The Lepton-specific 
model is similar to Type-I apart from the fact 
that the leptons couple to $i, instead of $ 2 . The 
Flipped model is similar to Type-II apart from the 
leptons coupling to 4> 2 , instead of $i. In all these 


models, the limit cos(/3 — a) —» 0 is such that the 
light CP-even Higgs boson, h, has indistinguish¬ 
able properties from a SM Higgs boson with the 
same mass. The cross sections for production by 
gluon fusion are calculated using SusHi [3iM7T| and 
the branching ratios are calculated with 2HDMC iooj . 
For the branching ratio calculations, it is assumed 
that tua = van = m#±, irih = 125 GeV and 
m\ 2 = m\ tan/3/(l + tan 2 /3). 

The constraints derived from the combined 
search in rr and bb final states are presented as 
a function of 2HDM parameters. The exclusion 
region in the cos(/3 — a) versus tan/3 plane for 
in a = 300 GeV are shown in Fig. [3] for the four 
2HDM models , while the constraints obtained in 
the m^-tan/3 plane for cos(/3 — a) = 0.10 are shown 
in Fig. [5] The width of the A boson in the 2HDM 
may be larger than the experimental mass resolu¬ 
tion, and it is taken into account in the 2HDM pa¬ 
rameter exclusion regions for widths up to 5% of 
m,A- For Type-II and Flipped models, Higgs boson 
production in association with 6-quarks dominates 
over gluon fusion for large tan /3 values (tan /3 > 10). 
The cross section for the 6-associated production 
uses an empirical matching of the cross sections in 
the four- and five-flavour schemes [58|. Cross sec¬ 
tions in the four-flavour scheme are calculated ac¬ 
cording to Refs. [S3 EH] and cross sections in the 
five-flavour scheme are calculated using SusHi. The 
relative efficiencies for the 6-associated and gluon 
fusion production as well as the predicted cross- 
section ratio are taken into account when deriving 
the constraints in the two-dimensional planes shown 
in Fig. [4] The 6-associated production efficiencies 
are estimated using PYTHIA8 and SHERPA samples. 
The regions of parameter space excluded at 95% 
CL by the A —>• rr decay mode are displayed in the 
same plots, using the results of a search for a heavy 
Higgs boson decaying into rr (Ref. |T3]), reinter¬ 
preted considering only the production of an A bo¬ 
son via gluon fusion and 6-associated production. 
For niA values below the ti kinematic threshold, 
the search presented here can exclude cos(/3 — a) 
values down to a few percent for tan /3 values up to 
« 3. 

7. Conclusions 

Data recorded in 2012 by the ATLAS experiment 
at the LHC, corresponding to an integrated lumi¬ 
nosity of 20.3 fb -1 of proton-proton collisions at 
a centre-of-mass energy 8 TeV, are used to search 
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for a CP-odd Higgs boson, A, decaying to Zh, 
where h denotes a light CP-even Higgs boson with 
a 125 GeV mass. No deviations from the SM back¬ 
ground predictions are observed in the three final 
states considered: Zh —>• ££tt, Zh —> Mbb , and 
Zh — > vvbb. Upper limits are set at the 95% confi¬ 
dence level for a{gg — > A) x BR(A —► Zh)x BR(ft, —> 
//) of 0.098-0.013 pb for / = r and 0.57-0.014 pb 
for f = b in the range of tua = 220-1000 GeV. This 
Zh resonance search improves significantly the pre¬ 
viously published constraints on CP-odd Higgs bo¬ 
son production in the low tan /? region of the 2HDM. 
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